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ABSTRACT 


neq  n  at  vapor  pressure 


Two  formulas,  one  for  the  characteristic  time  of 
filling  a  void  with  the  vapor  of  the  surrounding  liq¬ 
uid,  and  one  of  filling  the  void  by  diffusion  of  the 
dissolved  gas  in  the  liquid,  are  derived.  Based  on 
this  analysis,  it  is  seen  that  in  an  oil  film  bearing 
operating  under  dynamic  loads,  the  content  of  the 
cavitation  region  should  be  oil  vapor  rather  than  the 
air  liberated  from  solution,  if  the  oil  is  free  of 
entrained  air. 

NOMENCLATURE 

C  mass  concentration  of  gas 

Cs  saturation  value  of  C 

D  coefficient  of  binary  diffusion 

H  Henry's  constant 

h  thickness  of  liquid  film 

*cond  number  of  molecules  condensing  on  unit  surface 
area  in  uni t  t ime 


p  pressure 

ps  pressure  corresponding  to  Cs 

R  radius  of  circular  void 

r  radial  distance  in  plane  polar  coordinate  sys¬ 

tem,  measured  from  center  of  void 

,'Jt  universal  gas  Constant.  =  8.314  kJ/(kmole-k) 

T  absolute  temperature 

t  time  coordinate 

*diff  characteristic  time  of  diffusion,  Eq.  (22) 

*dyn  characteristic  time  of  dynamic  operation 

tgvap  characteristic  time  of  evaporation,  Eq .  (4) 

vm  mean  molecular  speed 

5  dimensionless  length  coordinate 


j d i f f  ®ass  of  gas  diffused  into  the  void  from  unit 

surface  area  of  liquid  in  unit  time 

j  number  of  molecules  evaporating  from  unit  sur- 

evap 

face  area  in  unit  time 


p  density 

t  dimensionless  time  coordinate 

Subscripts: 


k  Boltzmann  constant,  -  1 .380622x10-23  j /K 


g  of  gas 


8.  thickness  of  diffusion  layer 


8  of  liquid 


M  molecular  weight 


INTRODUCTION 


m  mass  of  a  molecule 

n  number  of  molecules  in  unit  volume 


The  cavitation  of  oil  films  in  bearings  operating 
under  dynamic  loads  is  not  well  understood.  One  of  the 
unresolved  questions  is  whether  the  cavitation  bubble 
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contains  air,  or  vapor,  or  is  just  an  empty  void.  This 
question  is  closely  tied  with  the  behavior  of  the  bub¬ 
ble,  e.g.,  whether  it  collapses  or  sustains  itself  dur¬ 
ing  the  high  pressure  portion  of  the  pressure  cycles. 
Understanding  the  bubble  behavior  is  indispensable  to 
the  construction  of  a  sound  theory  of  dynamically 
loaded  bearings.  At  present,  reliable  design  informa¬ 
tion  is  still  lacking  for  important  applications  of 
dynamically  loaded  bearings,  such  as  the  squeeze  film 
damper . 

Air  and  vapor  cavities  arise  in  oil  film  bearings 
in  several  ways  (Dowson  and  Taylor,  1979).  Air  from 
the  surrounding  atmosphere  may  be  entrained  in  oil  when 
subambient  pressure  occurs  in  the  bearing  clearance; 
the  dissolved  air  in  oil  may  be  liberated  when  the 
pressure  falls  below  the  saturation  pressure  of  the 
solution;  and  oil  molecules  may  evaporate  when  the 
pressure  falls  to  the  vapor  pressure  of  oil.  Air 
entrainment  is  device  dependent;  it  can  be  dealt  with 
only  according  to  specific  circumstances.  On  the  other 
hand,  liberation  of  dissolved  air  and  oil  evaporation 
are  processes  dependent  only  on  fluid  properties;  they 
can  be  treated  in  general  terms.  In  the  following,  the 
latter  two  processes  are  analyzed  and  their  respective 
time  scales  determined.  By  comparing  these  time  scales 
with  that  of  bearing  operation,  one  can  judge  the  rele¬ 
vance  of  these  processes  with  reference  to  the  cavita¬ 
tion  in  dynamically  loaded  bearings. 

CHARACTERISTIC  TIME  OF  EVAPORATION 

Consider  that  a  void  is  suddenly  created  in  a  liq¬ 
uid  film,  the  molecules  start  to  evaporate  from  the 
liquid  surface  into  the  void.  After  some  time,  the 
pressure  in  the  void  is  built  up  to  the  vapor  pressure 
of  the  liquid,  which  signifies  the  equilibrium  between 
the  flux  of  evaporation  and  flux  of  condensation.  The 
vapor  pressure  of  the  liquid  is  assumed  low.  Hence, 
one  may  consider  (Langmuir,  1913)  that  the  flux  of 
evaporation  is  independent  of  the  presence  of  vapor  in 
the  void:  and  that  the  flux  of  condensation  is  likewise 
not  affected  by  the  flux  of  evaporation,  and  is  deter¬ 
mined  only  by  the  pressure  in  the  void.  Since  the  flux 
of  molecules  is  one-fourth  the  number  density  times 
mean  molecular  speed  (Present,  1958.  p.  22),  one  may 
write: 

i  =  7  n  v  ,  j  ,  =  7  nv  (la.b) 

'evap  4  eq  m  'cond  4m  '  ’ 

Consider  that  the  film  has  thickness  h  and  the  void 
is  circular  with  radius  R.  Then,  the  change  in  the 
number  density  in  the  void  may  be  expressed  as: 


-t/t 

1  -  e  evap 


The  parameter  tevap  signifies  the  time  needed  to 
build  up  the  pressure  in  the  void  to  63  percent  of  the 
vapor  pressure,  and  may  be  taken  as  the  characteristic 
time  of  evaporation.  It  is  interesting  to  see  from 
Eq.  (4)  that  tevap  is  essentially  the  time  a  molecule 
takes  to  travel  across  the  void,  and  is  independent  of 
the  vapor  pressure  of  the  liquid.  With  the  mean  molec¬ 
ular  speed  given  as  (Present,  1958.  p.  78) 


which  is  based  on  the  Maxwell  distribution,  this  char¬ 
acteristic  time  can  be  readily  calculated. 

In  the  following,  an  oil-like  substance,  squa- 
lene,  of  known  molecular  structure  (Morales  and 
Buckley,  1988)  is  used  for  the  calculation.  The 
needed  properties  of  squalene  are  given  in  Table  1. 

From  the  observation  of  the  cavitation  in  a  dynami¬ 
cally  loaded  bearing  (Sun  and  Brewe.  1990),  the  cav¬ 
itation  region  was  large  and  covered  a  fraction  of  the 
bearing  surface.  If  the  typical  size  of  the  cavita¬ 
tion  bubble  is  taken  to  be  R  »  10  mm  and  the  operating 
temperature  to  be  T  =  293  K.  then  Eqs.  (6)  and  (4) 
yield  vn  =  120  m/s  and  tevap  =  0.167  ms. 

CHARACTERISTIC  TIME  OF  DIFFUSION 

The  solubility  of  gas  in  liquid  is  given  by 
Henry  s  law: 


C 

M 


which  states  that  the  mole  fraction  of  the  gas  in  the 
solution  is  proportional  to  the  pressure  of  the  gas 
outside  the  solution.  The  coefficient  appearing  in 
this  relation  is  called  Henry  s  constant  whose  value 
depends  on  the  particular  combination  of  gas  and  liq 
uid.  Henry's  law  applies  to  weak  solutions,  i.e., 


Hence.  Eq.  (7)  may  be  written  approximately  as: 


TfR2h  ~  =  2 trRh (  j  -  j  ,) 
dt  '’evap  ’cond’ 

Combining  Eqs.  (1)  and  (2),  one  obtains: 


where 


(2) 


(3) 


C 


Consider  a  liquid  film  saturated  with  a  dissolved  gas. 
According  to  Eq.  (9).  the  saturation  concentration  and 
the  corresponding  saturation  pressure  ps  are  related 
by: 


The  solution  to  Eq.  (3)  satisfying  the  initial  condi¬ 
tion  n (0 )  -  0  is: 


•> 


Let  a  circular  void  of  radius  R  be  suddenly  created 
in  the  film  and  the  dissolved  gas  start  to  diffuse  into 
the  void.  It  is  intended  to  find  how  the  pressure  in 
the  void  increases  with  time. 

The  variation  of  gas  concentration  in  the  solution 
is  governed  by  the  diffusion  equation: 


I  L  (.  1  3C  R 

r  3r  \  3ry  D  3t  '  K 


The  initial  condition  is: 


C(r.O)  -  Cs  (11) 

The  boundary  condition  at  large  distance  from  the  void 


C(=,t)  •*  Cs  (U 

With  pit)  denoting  the  pressure  in  the  void,  Eq.  (9) 
is  applied  as  the  boundary  condition  at  r  =  R.  let 
(dyn  be  a  representative  duration  for  the  occurrence 
of  cavitation,  and  4  the  characteristic  size  of  the 
region  where  diffusion  effectively  takes  place.  From 
Eq.  (10 1  one  deduces  that 


1  ■  ("‘dJ'" 


For  a  gas-in-liquid  system,  typically  D  -  0(10~9 
m-/s):  and  in  the  operation  of  dynamically  loaded  bear¬ 
ings,  typically  t^yn  -  0(10"2  s).  Then  Eq.  (13)  shows 
that  4  -  0(10~6  ra).  Hence,  diffusion  is  effective 
oniy  in  a  thin  layer  whose  thickness  is  much  smaller 
than  R.  Substituting 


r  •  R  +  4?,  5  ~  0(1) 


1  *  1  dynT ' 


T  ~  0(D 


(14a, b) 


into  Eq .  (10)  and  retaining  only  the  leading  order 
terms,  one  obtains: 

32C  3C  , 

a  =  a-t 

35“ 

The  initial  and  boundary  conditions  are  transformed 
into: 


At  x  -  0,  C  -  C„ 


As  5  ■*  C  *  C 


At  5-0, 


C_ 

JV.  .  B 


(16a, b,c) 


result,  which  will  show  that  p  is  indeed  a  slowly 
varying  function  during  the  diffusion  process. 

The  mass  flux  into  the  void  is  given  by  Fick  s 

law : 


■  D@) 


Using  Eq.  (17)  to  evaluate  the  concentration  gradient, 
one  finds: 

•di'f  ■  ^f1  -  jhfe)  119 

V  S/ 

The  change- of  the  gas  density  in  the  void  may  then  be 
expressed  as: 

irR  n  j“-=  2TrRhjdiff  120; 

It  is  assumed  that  the  gas  in  the  void  obeys  the  ideal 
gas  law: 


P  -  T 


In  addition,  a  parameter  tdif.f  is  defined  to  simplify 
the  notation: 


/  y  h  V  "> 

t  =  rU 

diff  "  161  o ,,^T  ;  D 

Substituting  Eqs .  (3a).  (131.  (21).  and  (22)  into  |20|, 
one  obtains  the  equation  governing  p: 

'  dB  _• _ -I-— _  |-3| 

dt  ,/  ',1/2  '-31 

diff1/ 

The  solution  to  this  equation  satisfying  the  initial 
condi tion  p(0)  =  0  i s : 


‘diff/ 


The  parameter  tjjff  signifies  the  time  needed 
for  the  pressure  in  the  void  to  reach  63  percent  of 
ps,  and  may  be  taken  as  the  characteristic  time  of  dif¬ 
fusion.  Its  evaluation  requires  the  knowledge  of  the 
coefficient  of  binary  diffusion  and  Henry  s  constant. 
For  air  in  squalene,  these  quantities  are  found  (Appen¬ 
dix)  to  be  D  -  0.495xi0~9  m2/s  and  H  =  13.4  MPa. 
Again  by  taking  R  «  10  mm  and  T  =■  293  K.  one  obtains 
from  Eq.  (22)  that  tdiff  =  0.324xl06  s. 

DISCUSSION 


The  solution  to  this  system  can  be  readily  found.  In 
terms  of  the  original  variables,  it  is: 


1  -  erf 


This  is  an  approximate  solution  because  during  its 
derivation  p(t)  is  treated  as  constant.  This  quasi¬ 
steadiness  approximation  can  be  justified  by  the  final 


The  characteristic  time  of  diffusion  calculated  in 
the  above  example  appears  unexpectedly  long,  and  war- 
rents  a  comparison  with  the  results  of  similar  studies 
in  the  literature.  Epstein  and  Plesset  (1950)  calcula¬ 
ted  the  time  of  completely  dissolving  a  spherical  air 
bubble  in  water.  In  their  analysis  the  concentration 
of  air  at  the  bubble  wall  was  assumed  always  saturated 
(Cs),  and  that  in  water  (C, )  waj  treated  as  a  parame¬ 
ter.  The  dissolving  time  ,*"’s  expressed  as  a  function 
of  the  concentration  ratio  (Cj/Cs)  and  the  initial 


uoe> 
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bubble  radius  ( R0 ) •  In  the  case  of  Cj  -  0,  where 
diffusion  was  the  fastest,  the  dissolving  time  was 
found  to  be  1.25  s  for  R0  -  0.01  mm;  1.25xl02  s  for 
R0  -  0.1  mm;  and  1.2Sxi04  s  for  R0  »  1  mm.  Recogni¬ 
zing  from  Eq.  (22)  that  tjjff  is  proportional  to  the 
square  of  the  bubble  radius,  one  would  expect  a  dis¬ 
solving  time  of  0(10®  s)  for  a  bubble  radius  of  10  mm. 
Krieger,  Mulholland,  and  Dickey  (1967)  actually  mea¬ 
sured  the  dissolving  time.  For  oxygen  bubbles  of  ini¬ 
tial  radii  in  the  neighborhood  of  0.3  mm,  they  observed 
the  dissolving  time  to  vary  from  6xl02  s  to  0.45xl04  s, 
for  Cj  ranging  between  zero  and  90  percent  of  satura¬ 
tion.  Thus,  the  credibility  of  Eq.  (22)  is  supported 
by  these  earlier  results;  they  all  point  to  the  fact 
that  diffusion  is  a  very  slow  process. 

The  load  frequencies  encountered  in  current  oil 
film  bearings  fall  in  a  certain  range.  The  operating 
speed  of  the  squeeze  film  dampers  in  gas  turbine 
engines  may  be  as  high  as  a  few  ten-thousand  revolu¬ 
tions  per  minute;  while  that  of  the  main  and  connect¬ 
ing  rod  bearing  in  reciprocating  internal  combustion 
engines  may  be  as  low  as  a  few  hundred  revolutions  per 
minute.  This  range  of  operating  speeds  corresponds  to 
a  t^yn  value  from  1  ms  to  0.1  s.  Thus,  tjyn  is  usu¬ 
ally  much  longer  than  the  characteristic  time  of  evapo¬ 
ration,  and  one  may  expect  the  presence  of  oil  vapor  in 
the  cavitation  bubble.  On  the  other  hand,  tjyjj  is 
much  too  short  compared  with  the  characteristic  time  of 
diffusion.  Hence,  there  is  not  enough  time  for  the 
dissolved  air  to  come  out  of  the  solution  and  find  its 
way  into  the  cavitation  bubble.  This  analysis  is  con¬ 
sistent  with  the  experimental  results  of  Sun  and  Brewe 
(1990),  where  no  difference  in  the  cavitation  pattern 
was  observed  between  the  degassed  and  the  aerated  oil 
throughout  the  cases  studied.  Since  the  diffusion  of 
dissolved  air  does  not  play  a  role  in  the  dynamic  cavi¬ 
tation  phenomenon,  then,  whether  the  oil  has  dissolved 
air  in  it  or  not  should  make  no  difference. 

It  has  often  been  reported1  that,  under  dynamic 
bearing  loads,  persistent  and  numerous  gas  bubbles  were 
found  in  the  oil  film,  which  eventually  became  a  well- 
mixed  two-phase  fluid.  The  presence  of  these  gas  bub¬ 
bles  has  been  attributed  to  both  the  liberation  of 
dissolved  gases  in  oil  and  the  entrainment  of  ambient 
air.  It  seems  clear  that  the  former  suspected  cause 
can  be  ruled  out  by  the  result  of  the  present  analysis, 

CLOSURE 


tour  to  the  Chemists'  wonderland  of  diffusion.  DC. 

Sun  is  particularly  grateful  to  Dr.  Edward  J.  Bissett 
of  the  General  Motors  Research  Laboratories  tor  his 
timely  objection  to  an  earlier,  incorrect  derivation 
in  the  work. 
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APPENDIX 

CALCULATING  D  AND  H  FOR  THE  SOLUTION  OF  AIR  IN 
SQUALENE 

Some  needed  properties  of  air  and  squalene  are 
listed  in  Table  1.  The  coefficient  of  binary  diffusion 
for  gas  in  liquid  may  be  calculated  by  using  the  Arnold 
equation  (Simons  and  Ponter,  1975)  which  reads: 


Two  formulas  (Eqs.  (4)  and  (22))  are  derived  which 
can  be  used  conveniently  to  determine,  respectively, 
the  characteristic  time  of  filling  a  void  with  the 
vapor  of  the  surrounding  liquid,  and  that  of  filling  a 
void  by  diffusion  of  the  dissolved  gas  in  the  liquid. 
These  time  scales  can  be  compared  with  the  characteris¬ 
tic  time  of  bearing  operation  to  determine  the  rele¬ 
vance  of  these  processes  with  reference  to  the  dynamic 
cavitation  in  oil  film  bearings. 

The  computation  shows  that  the  evaporation  process 
is  usually  fast  enough  to  fill  the  cavitation  bubble 
with  oil  vapor;  whereas  the  diffusion  process  is  much 
too  slow  for  the  dissolved  air  to  liberate  itself  and 
enter  the  cavitation  bubble. 
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where  the  empirical  constant  0.010  is  accurate  at 
20  °C;  the  units  of  the  quantities  on  the  RHS  of  the 
equation  are  as  shown  in  Table  1:  and  D  is  in  cm2/s. 
Substituting  the  appropriate  values  of  Table  1  into 
the  equation,  one  obtains  D  =  0.495xl0~°  cm2/s. 

Henry's  constant  can  be  estimated  by  the  method  of 
Osburn  and  Markovic  (1969).  The  procedure  is  to  first 
calculate  a  parameter  'parachor'  for  the  liquid  based 
on  its  molecular  structure.  From  the  parachor  and  the 
liquid's  molal  volume,  its  surface  tension  can  be  eval¬ 
uated.  Then  from  empirical  correlations  an  Ostwald 
coefficient  is  obtained,  which  can  be  converted  to 
Henry's  constant  in  unit  of  atmospheres.  Following 
this  procedure  and  with  the  use  of  Table  1,  one  obtains 
H  -  132  atmospheres. 
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TABLE  1.  -  SOME  PROPERTIES  OF  AIR  AND  SQL'ALENE 


Ai  r 

Squalene 

Molecular  weight 

Mass  of  a  molecule 
Density 

Molal  volume  (at  boil¬ 
ing  temperature) 
Absolute  viscosity 
Parachor 

Mg  =■  28.97  g/mole 

Vg  =  30  ml/mole 

M(  -  422 .8  g/mole 
m  -  70.2xl0-23  g 

Pt  =  0.81  g/ml 

Vt  =  522  mi /mole 

Pt  -  9-72  cP  (at  20  °C) 
[P]  -  1160.3 
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